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In the wake of a new kind of phase generally occurring in 
mesoscopic transport phenomena, we discuss the validity of 
Friedel sum rule in the presence of this phase. We find that 
the general Friedel sum rule may be violated. 

PACS numbers: 73.23.-b; 73. 23. Ad; 72.10.-d 



With large scale research in Mesoscopic Physics over 
the last few decades, many of the well established no- 
tions of Condensed Matter Physics has been found to be 
violated in mesoscopic samples. Breakdown of Onsager 
reciprocity relation , violation of Ohms law J 1 , absence 
of material specific quantities like resistivity |3 [, violation 
of Hund's rule Q etc., are a few such examples. The pur- 
pose of this work is to show the violation of Friedel sum 
rule in mesoscopic systems. 

Friedel sum rule relates the density of states inside a 
fixed potential scatterer to the scattering phase shifts || . 
A deduction of the sum rule can be found in many text 
books and intuitively it can be understood as fol- 
lows. Consider for example a fixed spherically symmetric 
potential scatterer. Now we enclose it in a larger spher- 
ical volume. In an energy interval dE, the number of 
states depend on the number of times the specific bound- 
ary conditions can be fulfilled by the wave function of the 
electron. So when the energy is changed, it introduces a 
phase shift of the electron wave function and so changes 
the number of times the specific boundary condition can 
be satisfied. And hence the density of states p' inside the 
impurity is related to the scattering phase shift r\ in the 
following fashion B. 



drj 
dE 



= 7T/9 



(1) 



This can be extended to the partial wave analysis of scat- 
tering states and many important issues can be under- 
stood in terms of the Friedel sum rule Q. In case of 
a non-spherical scatterer or non-spherical Fermi surface, 
the scattering matrix is in general an NxN matrix. For 
any such general NxN scattering matrix S, the Friedel 
sum rule can be written as [pfl 



dO/dE = irp', 



(2) 



where 8 = £^ £j, exp[2i^i] being the eigenvalues of the 
scattering matrix S. This can be further written in a 
compact form as 



LJL(ln(det[S}))=7rp'. 



(3) 



For one-dimensional systems where the scattering matrix 
is 2x2, the Friedel sum rule was thought to be further 
simplified to give 0] 



dargjt) 
dE 



7T/9 



(4) 



where t is the transmission amplitude but this is not true. 

Recently a new phase has been discussed in Ref. Q for 
scattering by a stub where the scattering matrix is 2x2, 
and it is believed that this phase is also observed 

in mesoscopic systems experimentally |D| . This phase 
is a general feature of transmission zeroes that always 
occur in Fano resonances in Quantum Wires and Dots, 
the stub structure being the simplest example fl2|,|l3| ■ It 
was shown in Ref. that the phase slips is a new phase 
associated with the violation of parity effect because it 
is different from Aharonov-Bohm phase, statistical phase 
and phase due to wave- like motion of electrons depending 
on their wave vector or energy. Had it not been different 
from the other three phases, parity effect would not have 
been violated [{Hi. The specialty of this phase is that it 
is discontinuous as a function of energy, i.e., the phase 
of the wavefunction changes by 7r although its energy 
does not change. To be more precise this phase does not 
originate from change in wave-vector due to change in 
energy. Hence in view of the discussions before Eq. 1 
one can question the validity of Friedel sum rule in the 
presence of this phase [^2|-pT|. We shall give a pictorial 
description of this phase later (short-dashed and long- 
dashed curves in Fig. 1). 

The scattering matrix for the stub is 



S = 



r t 
t r 



(5) 



where r and t are reflection and transmission amplitudes 
across the stub and are 



r = cos[kL]/ {—cos[kL\ + 2isin[kL]) 



(6) 



and 



t = {-2isin[kL]) / \cos[kL] - 2isin[kL]). (7) 
The eigenvalues of the S matrix are 

(cos[kL\ + 2isin[kL\) / {— cos[kL] + 2isin[kL]) and — 1 

(8) 

Hence as defined in Eq. 2 
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ArcTan[-4cos{kL}sin[kL]/(-cos[kL} 2 + Asin[kL} 2 )} 

(9) 
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In Fig. 1 we plot 9 (solid curve), argif) (short-dashed 
curve) and arg(r) (long-dashed curve) (given in Eqs. 
6, 7 and 9) versus kL. It can be seen that arg(t) 
and arg(r) show discontinuous jumps and drops by tt 
JiT| but they cancel in such a way that 9 is continuous 
and monotonously increasing. Hence one finds that the 
Friedel sum rule (Eqs. 2 and 3) is not violated al- 
though because of the discontinuous slips in arg(t) Eq. 4 
is obviously violated because density of states can never 
be infinite while the LHS of Eq. 4 can be infinite. And 
hence one can say that so far no one has found a violation 
of Friedel sum rule (Eqs. 2 and 3). We shall show the 
violation of the Friedel sum rule in the presence of this 
new phase. 

Transport across the stub structure has acquired a lot 
of importance recently p] , p^|JT4|JT^ | . All analysis so far 
are based on calculations with a hard wall boundary con- 
dition (an infinite step barrier potential or an infinite step 
well potential) at the dead end of the stub (we refer to it 
as the hard walled stub and for which Eqs. 6, 7, 8 and 9 
are derived). An infinite potential well at the dead end of 
the stub reflects an incident electron with unit probabil- 
ity. Now a small perturbation from this would be a finite 
but very deep potential well at the dead end of the stub 
(soft walled stub). Electrons are almost entirely reflected 
from the end of the stub and a negligible fraction escapes. 
Dephasing can also give similar escape probability. The 
scattering problem in this case is depicted in Fig. 2 and 
also explained in the figure caption. It is solved using the 
mode matching technique or Griffith's boundary condi- 
tions jl7]], that give the continuity of wavefunction and 
the conservation of currents at the junctions. In this case 
the transmission zero in x-direction is replaced by a mini- 
mum |l8[| . We first intend to understand what happens to 
the discontinuous phase change that occur due to trans- 
mission zeroes in this case. So in Fig. 3 we plot transmis- 
sion coefficient T = \t\ 2 (solid curve) and the argument 
of the transmission amplitude t (short-dashed curve) in x 
direction, versus kL for an almost hard walled stub. The 
transmission coefficient shows very deep minima and at 
the same points arg(t) show very sharp but continuous 
drops. For the completely hard walled stub there is an 
exact zero and associated with it a discontinuous slip 
by 7r as shown in Fig. 1. In the same figure (Fig. 3) 
we also plot transmission coefficient in the x-direction 
(dash-dotted curve) and the corresponding argument of 
the transmission amplitude (long-dashed curve) versus 
kL for a very soft walled stub. At the points where the 
solid curve show very deep minima, dash-dotted curve 
show shallow minima. Also the fast phase drops change 
over to a slower decrease. 

Having understood the phase slips further we move on 
to the three prong scatterer (Fig. 4) that is often encoun- 
tered in mesoscopic systems [[L9| including the experimen- 
tal set up of Ref. ]l5[ | and many such similar experiments. 
The scattering problem in this case is described in the fig- 
ure caption. From the continuity of wavefunctions (first 
Griffith's boundary condition) we get the following equa- 



tions (variables and parameters are defined in Fig. 4 and 
it's caption). 

1 + r = a exp[—iqLi] + b exp[iqLi] ; a + b = c + d: 

a + b = f + g, cexp[iqL2\ + d exp\—iqL-i\ = e; 

/ exp[iqL 3 ] + g exp[-iqL 3 ] = h. (10) 

And from the second Griffith's boundary condition which 
is the conservation of currents at the junctions (S^^ 1 = 
0, that can be derived from current conservation, here ipi 
is a wavefunction at a junction, Xi is coordinate at that 
junction, and the sum over i stands for all such wavefunc- 
tions incoming or outgoing at a junction, the convention 
followed is that currents flowing into the junction is posi- 
tive while currents flowing out of the junction is negative) 
we get the following equations. 

k — kr — qa exp[—iqLi] + qb exp[iqLi] = 0; 

a — b — c + d — f + g = 0; 

qcexp[iqL2\ — qdexp\—iqL-^ — fee = 0; 

qf exp[iqL 3 ] - qg exp[-iqL 3 ] -kh = 0. (11) 

Thus we have 9 equations and exactly 9 unknown quanti- 
ties (a, 6, c, d, e, /, g, h and r) and so the problem is com- 
pletely defined. Once the unknowns are solved, the wave- 
function is known at all points exactly and so the density 
of states as well as the scattering matrix can be calcu- 
lated exactly. The scattering matrix in this case is 

(tn ti2 ti 3 \ 
fai h2 t 23 J . (12) 
*31 ^32 ^33 / 

Here in=r=transmission amplitude to the first prong 
when the incident beam is from the first prong. ti% = e 
is the transmission amplitude to the second prong when 
the incident beam is from the first prong. £13 = h is 
the transmission amplitude to the third prong when the 
incident beam is from the first prong. The other matrix 
elements are to be calculated when similar incident beam 
of unit flux is from the other two directions in Fig. 4. For 
the case of Fig. 4, the partial density of states is given 
by the following expression 

pi = irpi = - — [ / \aexp[iqx] + bexp\—iqx]\ 2 dx+ 
hv J- Ll 

rL-2 

I \cexp[iqy] + dexp[— iqy] \ 2 dy+ 
Ja 

/ \fexp[iqz] + gexp [-iqz]\ 2 dz], (13) 
Jo 
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where v = hk/m. a, 6, c, d, f and g are determined from 
Eqs. 10 and 11. P2 and pz are to be evaluated similarly 
when the incident beam is from the other two directions 
in Fig. 4, and p = pi + P2 + P3- For the symmetric three 
prong scatterer (L± = L2 = L3 = L), the antircsonances 
are almost cancelled by the resonances but still a vio- 
lation of Eq. 2 or 3 can be seen at low energy. This 
cancelling effect of resonance and antiresonance can be 
avoided by choosing incommensurate values of {L\ + L3) 
and L2, i.e., for asymmetric configurations. We will now 
go to the asymmetric configuration and demonstrate a 
large difference between r = 39 /dE — ^-J^[ln[Det[S]] 
and p = irp' at large energies (E « V). This is shown 
in Fig. 5. We want to emphasize that at very high en- 
ergy, compared to the energy scale V in the system, when 
multiple scattering and the new phase becomes insignif- 
icant, we recover Friedel sum rule perfectly. But when 
this new phase is present at energies (E < V), there is 
a large difference between r and p and hence a complete 
violation of Friedel sum rule. In Fig. 5, r or d9/dE, can 
become substantially negative, i.e., 9 can undergo a drop 
like arg(t) in Fig. 3. The new phase need not always 
appear as a drop but can also appear as a rise and then 
the LHS of Eq. 3 can remain positive all the time while 
deviating from the RHS of Eq. 3. This is shown in Fig 
6. 

Thus our exact calculation of density of states 
and scattering matrix elements show the deviation of 
^rJ^ln[Det[s\\ from 717/ in the presence of phase slips. 
The phase slips are a general feature of Quantum wires 
with defects jL2| and Quantum Dots and these phase slips 
are at the origin of drops in 9 and hence deviation or vio- 
lation of Friedel sum rule. Only 2x2 S-matrix is a special 
case where as shown in Fig. 1 some scattering matrix el- 
ements undergo a phase jump and some undergo a phase 
drop in such a manner that they cancel and the phase 
slips do not affect ln[Det[S}] or 9. The general feature 
is that they do not cancel and Friedel sum rule gets vi- 
olated. The attractive potential in the three prong scat- 
terer offsets the symmetry between the phase jumps and 
the phase drops so that they do not cancel each other. 

The author acknowledges useful discussions with Prof. 
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Figure captions 
Fig. 1 Arg(r) (long-dashed curve), arg(t) (short-dashed 
curve) and 9 (solid curve) for the hard walled stub. 
Length of the stub is L and it is taken to be the unit 
of length. We choose h = 2m = 1 . 

Fig. 2 A scattering problem with conventional notations 
is depicted here, k = \fE is the wave vector in the thin 
regions where the Quantum Mechanical potential is 0. 
q = s/E + V is the wave vector in the thick regions where 
the Quantum Mechanical potential is — V. x and y are 
coordinates and the origin of coordinates is also depicted 
in the figure, t and c are transmission amplitudes in x 
and y directions, respectively, while r is the reflection 
amplitude. Distance between points P and Q is L. 
Fig. 3 The solid curve is transmission coefficient T = \t\ 2 
across the soft walled stub described in Fig. 2. The short- 
dashed curve is the phase of the transmission amplitude 
t across the stub. We choose VI? = — 10 6 so that it 
is in the hard wall limit, and Ti = 2m = 1. Next we 
make VL 2 —-100 and plot the transmission coefficient T 
in dash-dotted curve. The phase of the transmission am- 
plitude t is given by long-dashed curve. 
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Fig. 4 A scattering problem with conventional notations 
is depicted here, k — \[E is the wave vector in the thin 
regions where the Quantum Mechanical potential is 0. 
q = yjE + V is the wave vector in the thick regions where 
the Quantum Mechanical potential is —V. x,y,z,u,v and 
w are coordinates and the origin of coordinates is also 
depicted in the figure, e and h are transmission ampli- 
tudes in respective directions, while r is the reflection 
amplitude. Distance between (u=0) and (x=0,y=0,z=0) 
is L\. Distance between (v=0) and (x=0,y=0,z=0) is L 2 . 
Distance between (w=0) and (x=0,y=0,z=0) is L 3 . 
Fig. 5 p (solid curve) and r = |§ =LHS of Eq. 3 (dotted 
curve) versus kL for the scattering problem described in 
Fig. 4. We choose VL 2 = -100, L x = L 3 = L, L 2 =AL 
and h = 2m = 1. 

Fig. 6 p (solid curve) and t = || =LHS of Eq. 3 (dotted 
curve) versus kL for the scattering problem described in 
Fig. 4. We choose VL 2 = -100, L 1 =L 3 = L, L 2 =2AL 
and h — 2m = 1. 
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